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Theory of the birefringence of the refractive index in atomic diamagnetic dilute gases in the presence of 
static electric (optical Kerr effect) and magnetic (Cotton-Mouton effect) fields is formulated. Quantum- 
statistical expressions for the second Kerr and Cotton-Mouton virial coefficients, valid both in the low 
and high temperature regimes, are derived. It is shown that both virial coefficients can rigorously be 
related to the difference of the fourth derivatives of the thermodynamic (pressure) virial coefficient with 
respect to the strength of the non-resonant optical fields with parallel and perpendicular polarizations 
and with respect to the external static (electric or magnetic) field. Semiclassical expansions of the Kerr 
and Cotton-Mouton coefficients are also considered, and quantum corrections up to and including the 
second order are derived. Calculations of the second Kerr and Cotton-Mouton virial coefficients of the 
■^Ho gas at various temperatures are reported. The role of the quantum-mechanical effects and the 
convergence properties of the semiclassical expansions arc discussed. Theoretical results are compared 
with the available experimental data. 



Keywords: Optical Kerr effect, Cotton-Mouton effect, quantum-statistical theory, semiclassical 
expansion, collision-induced properties, helium gas 



1. Introduction 

Atomic gases composed of closed-shell diamagnetic atoms are optically isotropic. This 
means that the speed of light traversing the gas sample is independent of the polarization of 
the light. External fields, like the electric or the magnetic fields, strongly modify refractive 
properties of gases leading to the optical anisotropy of the gas in the field. The optical 
birefringence of the refractive index of isotropic gases in the electric field was first observed 
by John Kerr in 1875 [1]. Kerr discovered that the refractive coefficient of a gas changes 
depending on the polarization of light with respect to the direction of the external static 
electric field. He found that the difference between the refractive coefficients for the light 
with parallel and perpendicular polarizations with respect to the field vector, n\\ and n_|_, 
respectively, is proportional to the square of the electric field £: 

n\\-n^ = Ke^, (1) 



where the proportionality constant has been referred to as the Kerr constant. In 1955 
Buckingham [2, 3] slightly modified the expression for the Kerr constant, essentially by 
taking the limit of the above expression to the zero field. We will come back to this point 
in sees. 2.1 and 4.4. 

First experimental observations of the birefringence of the refractive properties in the 
magnetic field were reported in 1907 by Aime Cotton and Henri Mouton [4]. The results 
of the experimental measurements show that similarly to the Kerr effect the difference 
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between the refractive coefficients for the hght with parallel and perpendicular polarization 
with respect to the magnetic field vector is proportional to the square of the applied field 
B: 

ny -n^ = CB"^, (2) 

where the proportionality constant has been referred to as the Cotton- Mouton constant. In 
1956 Buckingham and Pople [5] modified the expression for the Cotton-Mouton constant 
in the spirit of the modified definition of the Kerr constant, i.e. essentially by taking the 
limit of the above expression to the zero magnetic field. 

It was observed experimentally that at very low gas number densities p the Kerr constant 
Km depends linearly on p, and the proportionality coefficient is proportional to the atomic 
second hyperpolarizability. At higher pressures departure from the ideal gas law is observed, 
and it was shown that terms quadratic, cubic, and higher in p contribute to [2, 3]. 
For dilute gases the quadratic term is dominant, and is referred to as the Kerr virial 
coefficient. Buckingham and collaborators were the ffi'st to study, both theoretically and 
experimentally, the pressure effects on the Kerr constant [2, 3, 6, 7], and derived [2, 3] a 
classical expression for the Kerr virial coefficient in terms of the interatomic interaction 
potential and interaction-induced electric properties [8]. Identical behavior as function of 
p was observed for the Cotton-Mouton constant Cm, and the classical expression for the 
Cotton-Mouton virial coefficient was derived in Ref. [5]. 

It is well known that at very low temperatures thermodynamic and dielectric properties of 
gases depart from the classical picture. This was demonstrated in the classical works on the 
pressure virial coefficient of the helium gas at low temperatures [9-11] and in Refs. [12-15] 
for the virial expansion of the dielectric Clausius-Mossotti function. In particular Ref. [12] 
contains a very detailed discussion of the quantum effects on the dielectric virial coefficient, 
while Ref. [14] reports the semiclassical expansion of this coefficient to the second order, i.e. 
including the effects of the order of fi^. It was shown [12] that for temperatures above 100 
K the classical and quantum results differ by 2% at most for the helium-4 gas. At lower 
temperatures this deviation becomes larger and larger, and the semiclassical expansion 
diverges. 

Surprisingly enough, the refractive properties of gases were not studied with a quantum- 
statistical approach. Bruch and collaborators [13, 16] gave a quantum-statistical expression 
for the upper bound to the Kerr virial coefficient. Actually, the expression reported in Refs. 
[13, 16] cannot be correct since it contains singular objects like the squares of the quantum- 
mechanical operators. Rizzo and collaborators [15] noticed that the classical expressions 
for the dielectric virial and Kerr virial coefficients are very similar. In fact, one can obtain 
the expression for the Kerr virial coefficient by replacing the trace of the collision-induced 
polar izability tensor in the equation for the dielectric virial coefficient by a proper linear 
combination of the square of the collision-induced polar izability anisotropy and collision- 
induced trace of the second hyperpolarizability. Using this observation the Authors of 
Ref. [15] suggested that a proper quantum-statistical formulation of the Kerr effect can 
be obtained from the quantum-statistical description of the Clausius-Mossotti function 
[12] by simply replacing the collision-induced properties like in the classical case, but no 
mathematical proof that this is a correct procedure was given. In fact, in the present 
paper we show that the ad hoc procedure adopted by Rizzo and collaborators [15] is not 
correct. The very same remarks as above apply to the quantum-statistical treatment of 
the Cotton-Mouton effect reported in Ref. [15]. 

Given the fact that no systematic quantum-statistical studies of the refractive properties 
of the atomic gases in the external electric or magnetic fields are available in the literature, 
and that virial expansions at low temperatures are also useful in the field of the gas-phase 
NMR spectroscopy, as illustrated among others in Ref. [17] on the example of Arll2 for 
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which accurate theoretical data are available [18-20] , in the present paper we fill this gap 
and report systematic derivations of the quantum-statistical expressions for the Kerr virial 
and Cotton-Mouton virial coefficients. Our derived formulas are valid both in the low and 
high temperature regime. We also derive expressions for the quantum corrections in the 
semiclassical expansion up to and including the second-order h'^ term. The plan of this 
paper is as follows. In sec. 2 we report an expression connecting the Kerr virial coefficient 
to the fourth derivative of the thermodynamic (pressure) virial coefficient in the combined 
non-resonant and static fields. By using the quantum expression for the thermodynamic 
virial coefficient we derive an equation for the Kerr virial coefficient in terms of the collision- 
induced anisotropy of the polar izability tensor, collision-induced second hyperpolarizability, 
and eigenvalues and eigenfunctions of the Hamiltonian describing relative nuclear motion 
of two atoms in an interatomic potential. In this section we also present a systematic 
derivation of the semiclassical expansion of the Kerr virial coefficient, and report formulas 
for the first and second quantum corrections to the pure classical result. Sec. 3 is devoted to 
the Cotton-Mouton effect. In this section we briefly sketch the derivation of the quantum- 
statistical expression for the Cotton-Mouton virial coefficient, and report the final formula 
in terms of the electric and magnetic collision-induced properties. In sec. 4 we report 
numerical results illustrating our theoretical findings. All calculations will be reported 
for the helium-4 gas which shows the most pronounced quantum behavior in the low 
temperature regime. Our calculations will be based on the most recent ab initio potential 
for the helium dimer [21], on the anisotropy of the collision- induced polarizability tensor 
of Ref. [14], and on the collision-induced second hyperpolarizability and collision- induced 
magnetic properties of Refs. [22, 23]. It was argued in Ref. [24] that the 1996 results for 
the anisotropy of the collision-induced polarizability tensor [14] are not as accurate as 
those reported in Refs. [22, 24]. However, the data of Ref. [14] were shown to perfectly 
reproduce very precise measurements of the polarized and depolarized collision-induced 
Raman spectra in the high and low temperature regimes [25-27]. No proof of accuracy 
of the anisotropy models of Refs. [22, 24] by comparison with the experimental data was 
reported in the literature thus far. Therefore, in this paper we adopted the 1996 model of 
the collision-induced anisotropy of Ref. [14]. Finally, sec. 5 concludes our paper. 



2. Quantum-statistical theory of the optical Kerr effect 
2.1. Introductory remarks and definitions 

We consider a dilute gas composed of diamagnetic closed-shell atoms in a static electric 
field £ directed along the Z axis of the space-fixed coordinate system. In the presence of a 
non-resonant field with polarization parallel or perpendicular to the static field £, denoted 
and F"*-, respectively, the system shows a birefringence of the refractive coefficient. The 
Kerr coefficient Km is a direct measure of this optical birefringence in the limit of a small 
static field £. According to Buckingham [2, 3] the Kerr constant Km characterizing the 
optical birefringence of the refractive coefficient is defined by the following expression: 

Km = , , t ^ (3) 

(n2 + 2)2(e + 2)2^,->o £f ' ^ ^ 

where n and e are the refractive and dielectric constants of the gas in the absence of any 
external fields, respectively, ny and n± are the refractive constants of the gas in the presence 

of the non-resonant electric fields and F-^, and £i is the local field acting directly on the 
atoms in the gas, and thus different from the external static field £. By inserting into Eq. 
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(3) the Lorentz equation for the local field valid for a gas composed of spherical particles: 

£i = , (4) 

we arrive at the following expression for the Kerr constant: 

2n n\\ — re_L 

^-=3RWi-o^^- 

We will assume that the applied fields are small enough, so that ny + n_L ~ 2n. In such a 
case the expression for the Kerr coefficient, Eq. (5), becomes: 

i^. = -hm . (6) 

This expression was used in the past by Buckingham [2, 3] to derive the virial expansion of 
the Kerr constant in the high temperature regime, and will also be used in all derivations 
reported in the present paper. 



2.2. Virial expansion of the Kerr coefficient 

The Kerr constant is a function of the temperature T and gas number density p. We will 
now show that it can be represented as a virial expansion in the powers of p: 

Km = Akp + Bk{T)p^ + Ck{T)p^ + • • • , (7) 

where Ak is an atomic term independent of the temperature, and Bk(T) and Ck{T) are 
the second and third Kerr virial coefficients, respectively. We will prove that the power 
series expansion of the Kerr coefficient (7) is indeed correct and does not include any other, 
for instance, logarithmic or inverse power, dependence on p. To this end we recall the reader 
that for diamagnetic gases the refractive and dielectric constants are not independent, but 
are related by the following expression: 

V? = e. (8) 
We also know the virial expansion of the Clausius-Mossotti function [12]: 

^-^=A,p + B,{T)p^ + a{T)p^ + ■■■ , (9) 

where A^^ is related to the atomic polarizability oq by the expression: 

A. = (10) 

and an explicit expression for Bf^{T) valid in any temperature regime is known [12] in terms 
of the thermodynamic (pressure) virial coefficient in a static electric field £: 
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Here, E denotes the static electric field, is the Boltzmann constant, and B2{T;£) is the 
thermodynamic (pressure) virial coefficient at a temperature T and field £. 

By virtue of Eq. (8), we can rewrite the numerators in Eq. (6) in terms of the parallel 
and perpendicular Clausius-Mossotti functions: 

"1-1 .,1-1 ni-l a-1 



n 



2 



en +2' n2 +2 e_L + 2' 



where ey and ej_ are the dielectric constants of the gas in the presence of the external 
non-resonant fields, and F-^, respectively. By inserting the virial expansions of the 
Clausius-Mossotti function, Eq. (9), in the non-resonant fields parallel and perpendicular 
to the static field £, and F-*-, into Eq. (5) combined with Eq. (12): 

'1^ = 4{£)p + b!{T;£)p' + ..., (13) 
eil + 2 



e±-l 
e± + 2 



(14) 



and taking the limit <5 — we arrive at Eq. (7) with the Ak and Bk{T) coefficients defined 
as: 



A 



K 



18 



d^Al 



£=0 



d^Aj 



£=0 



(15) 



B 



K 



2TrkBT 
2^ 



d^B2{T-F\\-£) d^B2{T-F^-£) 



d£^dF^ 



d£^dF'i 



e=F«=F^=0 



(16) 



To derive the above equations we have used Eq. (11) and the fact that the dielectric 
constant e is an even function of the static external field £ [12], so it must be for the 

coefficients aI^^{£) and s|/^(r,f ): 



4/^{£) = AM + 



1 



d£'^ 



£=0 



£^ + 0{£' 



(17) 



i?|/^(r,£:) = s,(r,o) + 



IIA 



d£'^ 



£■=0 



£'^ + 0{£^). 



(18) 



Obviously, Ak can be written in terms of the atomic second hyperpolarizability 70: 



Ak = ^70. 



(19) 
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2.3. Quantum- statistical expression for the Kerr virial coefficient 

We start the derivation with the quantum-statistical expression for the thermodynamic 
(pressure) virial coefficient [28]: 

B2{T- F; £) = -\l dIl[W{K; F; £)-l], (20) 

where (R; F; £) is the Slater sum [28] in the presence of two fields, non-resonant F, and 
static £: 

W{n;F;£) = A|(R|e-^^(^'^)|R) + ^^-^{-R\e-^^(F'^)\R) , (21) 

j3 = l/ksT, Xb = {A-Kh'^/kBTm)^/^ is the thermal de Broglie wave length, / is the nuclear 
spin, and m the atomic mass. In view of the relation (20), Eq. (16) can be rewritten as: 



BKiT) = ^ /dR 



9^P^(R;F||;g) d^W(R;Fr,£) 
d£'^dF^ d£'^dF'j 



(22) 



Assuming only two-body interactions between the atoms in the gas, the Hamiltonian in 
the presence of two fields can conveniently be written as: 

#(F||/^ ,£) = + ^(1) (F||/^ , £) , (23) 

where i?^*^^ describes the relative motion of the nuclei in the absence of the fields: 

^(0) = _:?^V?> +y(i?), (24) 
m 

and V{R) is the interatomic interaction potential in the absence of any external fields. The 
field-dependent term H^^') is given by: 

H^'H^,£) = -lafAn)£' - la%{K)Fj - \l¥zM^)£'Pv (^5) 

and 

H'-'\F^,£) = -la%^zm£' - l<^¥xi^)Fl " ^7ll,xx(R)^^i^l- (26) 

The space-fixed components of the collision-induced polarizability tensor a^z and a^x 
can conveniently be written in terms of the body-fixed components and spherical angles 
{9, (p) of the R vector: 

afziK) = a{R) + ^ Aa(i?)P2(cos 6), (27) 



af^(R) = a{R) + ^a{R) | -L [C2,_2(e, <A) + C2,2(^,0)] - ^^'2(cos0)| 



(28) 



where Ci^m{9,(p) denotes the spherical harmonics in the Racah normalization, P;(cos^) is 
the Legendre polynomial, while a{R) and Aa{R) are the trace and anisotropy, respectively. 
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of the collision-induced polar izability tensor in the body-fixed frame: 

a{R) = ^[a^!{R) + 2a^!{R)], (29) 



Aa{R) = a^!{R)-a^!{R). (30) 

We assume here that the body-fixed z axis lies along the molecular axis. The components 
of the collision-induced polarizability tensor appearing in the expression above are defined 
as [8]: 

where is the component of the polarizability tensor of the dimer AB, while af^ and 
afj are the components of the polarizability tensor of the monomers A and B, respectively, 
all in the body-fixed frame. Similarly, the space-fixed components of the collision-induced 
hyperpolarizability tensor and zz conveniently be written in terms of the 

body-fixed components and the spherical angles (9, (p) of the R vector: 

+^ [^iZzziR) - ^^iZM + 247.^i^..(^)] C4,o(e, 0), (32) 



7lz,xx(R) 



1 

15' 



,BF 



21 



bzlAR) + ixlAR) - ^ixl-xiR)] C2A0, 



+ ^ bfF,zziR) + Ixl^AR) - ^^Ixl^xiR)] <P) + C2,-2(^, ' 

^ HiW,.AR) + ^^ixlzziR) - ^'^ixl^xiR)] CMO, <t>), 



where the three independent components are given by: 



7x 



^yy^yy — 37a;a;,' 



TO' 



Ix 



'~1yyy 



7. 



zz,zz^ 



(33) 



(34) 



and the single invariant of the second hyperpolarizability tensor can be written in terms 
of the body-fixed components as follows: 



7(^) 



^ zz,zz^R^ ^l'xx,yyiR) ~^ '^'yxx,zziR) 



„BF 



BF 



(35) 
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Figure 1. Integration contour in Eq. (37). 



Note that the above expressions are strictly valid for molecules of the D^oh symmetry. 

Our goal is to derive an expression for Bk is terms of the eigenvalues and eigenfunctions 
of fl'(o). To this end one has to perform the differentiation in Eq. (22). However, since 
the operators H^^'^ and H^^^ do not commute, the standard expression for the fourth 
derivative of the exponential function does not apply in this case. To simplify the notation 
let us rewrite H^^^ in the following symbolic form: 



2 2 4' 



(36) 



where the meaning of the operators A, B, and C is obvious from Eqs. (25) or (26). To 
derive an expression for the derivative of the exponential operator exp{—/3H) we make use 
of the following integral representation [29] : 



-fiH 



1 



dE- 



-I3E 



'2Tri J E-H 



(37) 



where the integration is done over the contour presented in Fig. 1. We can now expand 
the denominator appearing in Eq. (37) in the following way: 



1 



1 



E-H E-H(o)-HW {E-Hio)){l-{E-Hio)yim)) 

= (^-^(°))(T- (£;-#(°))-i^W)j"^ 

= (l-{E- ^W)-i#(i))-i(i? - #W)-i 



(38) 



1 



1 



l-{E-H(0)y^Hi^) E-H(^) 



Ef— 



n=0 



E - if(0) 



E - Hio) ' 
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and rewrite Eq. (37) as: 



1 /■ °° 1 



27ri / 



n=0 



1 



E - Hi^) 



(39) 



Note that the field dependence on the r.h.s. of the above expression appears solely in 
H^^\ so the exponential became a polynomial in the field strengths as variables, and the 
differentiations with respect to F and £ can easily be done: 



dE- 



-I3E 



-c- 



1 



e=F=o 2TTi J E - E - 



+ P{A, B)— i dE 



-B- 



1 



(40) 



2m J E-H(^) E-Hio) E - H^^) 



where P{A^ B) denotes the permutation of the symbols A and B. Assuming that we know 
the complete set of the eigenstates of H^'^^: 



H^^Mn)=E„ 



with the following wave function representation: 



n) 



^'„(R) = (R\n), 



and the following resolution of identity: 



(41) 



(42) 



(43) 



the expression for Bk{T), Eq. (22), becomes: 

Bk{T) ='^^e-'''-('iA-B-WJ...m - 7f;,.(R))l*»(R)> 



E 



En — Ek 



(^„(R)|aff(R)|*fc(R)) 



X (^fc(R)l(aff (R) - aff (R))|^„(R)) 



+ 



(44) 



27 21 + 



- e-'^^"(M/„(-R)|(7f,%(R) - 7.T..(R))l*n(R)) 
(v&„(-R)|aff(R)|^,(R)) 



n,k 



En — El 



X (^,(R)|(aff (R) - a^,^ (R))|*.(R)) 



SF/ 



When deriving the above expression we have used the following integral identities that 
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easily follow from the residue theorem: 



AEj^-^ = 27ri(-/3e"^^"), (45) 



dE— ^-^tttt; ^ = 2TTi 



{E-EnY{E-Ek) V {Ek-Er,Y En-Ek 



(46) 



At this stage we can use in Eq. (44) the representations (27)~(28) and (32)-(33) of the 
space-fixed components of the polarizability and hyperpolarizability tensors, respectively, 
and the following explicit forms of the bound and continuum eigenfunctions of H^'^^: 

^nJM(R) = ^^[^lj,M(e,.^), (47) 



^fcJM(R) = "i^^Yj^MiOA). (48) 
where ip„j{R) and ip^jiR) are solutions of the following radial Schrodinger equation: 

+ (y{R) + -^-EH = ^^ (49) 



m dB? \ m B? 

subject to the following normalization conditions: 



poo 

/ dRrnj{R)i^n,j{R) = 5nn,, (50) 

Jo 



POD 

/ dR rkAR)i'k,.j{R) = 6{Ek - Ek,). (51) 
Jo 

It is convenient to split the expression (44) into terms related to Aa and 7: 

Bk(T) = BkMT) + Bk^T), (52) 

where BK,-y{T) is given by the first and the third terms of Eq. (44), and Bx,Aa(T) by the 
remaining ones. 

Let us first consider the term Bx^-yiT). By inserting the explicit expressions for ^'„ja/(R) 
and ^fcjj\/(R), Eqs. (47)-(48), and using the following relation, 

*nJM(-R) = (-1)^*„JA/(R) = (-l)-^^^^y7,A/(e, 0), (53) 
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we arrive at the following expression for Bx .y{T): 

B,„(r) = ^i:(2j + i)(i + i^j 



J=0 



X 5]e-^^-^(Vnj(i?)|7(^)IV'nj(/2)> 

'- n 

poo 

Jo 

where the matrix elements appearing in the expression above are given by: 

POO 



(54) 



(55) 



The derivation of an explicit expression for the second term BK,Aa{T) is somewhat more 
involved and requires some angular momentum algebra. The final result reads: 



BK,Aa{T) 



2vrA| °° 



■n,n' 



J=0 J'=J,J±2 
^-l3E„j _ g 



En' J' — En J 



y y ( 1 + Lil:^:^ ) (2 J + 1) (2 J' + 1) ( " ^ 

405 4-: . 4^ , A 21 + 1 J ^ ^ 

-\{i;nj{R)\^a{R)\iJn'J'{R))\^ 



J J' 2 



(56) 



+ 2 J]/ dE, 



Ek — E^ 



nj 



\{i;E,j'iR)\Aa{R)\i;njiR))? 



dEk / dEk- 



Ek - Ek' 



\{^Ea{R)\Aa{R)\^PE,,J,iR))\' 



where the expression in the curly brackets is the 3j symbol [30]. The final quantum- 
statistical expression for Bk{T) is given by the sum of Eqs. (54) and (56). 



2.4. Semiclassical expansion 

The quantum-statistical expression for the Kerr virial coefficient is quite complicated. How- 
ever, at reasonably low temperatures, e.g. the liquid nitrogen temperature, the semiclassical 
expansion in powers of h"^ should work. To derive the semiclassical expansion of the Kerr 
virial coefficient let us recall the expansion of the Slater sum in powers of h? [31]: 

The dependence on the electric fields fI'/-*- and £ in the above expression can only enter 
through the field dependence of the interatomic potential V. The second derivative of the 
interatomic potential with respect to the field gives the interaction-induced polarizability, 
while the fourth derivative the interaction-induced second hyper polarizability [8, 15]: 

f d'V{R;Fr,£) \ ( d'V{R;F\\;£) \ 

[ dFl J, , ( dF^ ="^^^''^' ^''^ 

^ ^ F±=£=0 \ II / F||=£-=0 



11 



Thus, by inserting the expansion (57) into Eq. (22) we obtain the semiclassical expansion 
of the Kerr virial coefficient. It is convenient to keep the spUtting of Bk{T) into parts 
related to 7 and Aa: 



<Aa = <^A. + <]Aa + ^^A^- (61) 

The classical term of the zeroth-order in h? is given by the following expression derived 
by Buckingham in 1955 [2, 3]: 

, , o_2 foo 

<7(^) = |^/ dRl{R)e~^^^^^R', (62) 



B 



(0) 

K,Aa 



(T) 



8^ 
405 



dR (Aa(i?))'/3e-^^Mi?2. 



(63) 



The expression for the first quantum corrections B^j^^^{T) and -Sj^^^„(r) are somewhat 
more complex and involve the first derivatives of the potential and collision-induced prop- 
erties with respect to R: 



2 ^2/32^2 



243 



m 



dR 



(3- 



^dV d7 
^dRdR 



-PV(r)j^2^ (64) 



BKAc.iT)--———j^ dR- 



dAa 
dR 



Aa{R) 
R 



4Aa- 



dAady 



dR dR 



/3 



(2) 

The expression for the second quantum correction related to 7, Bj^^^{T)^ reads: 
,(2) 2 h^f3^7r' 



B^m 



243 m2 



di?[7(fi)/i(fi)/3 + /2(i^)]e-^^('■)i^^ 



(65) 



(66) 



where the auxiliary functions fi{R) and f2{R) are given by: 



/i(K) 



5R^\dRj 9R\dRj 72\dRj 



■1 ^-i^i^] + 



10 Vdi?2 



(67) 
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•'^^ ' 5R-^dRdR 3RdR\dRj 18dR\dRj ^ 

1 d^-f dV 
~ 5d¥dR^' 



(68) 



Finally, the contribution to the second quantum correction related to the anisotropy of the 
collision-induced polar izability is given by: 



B 



(2) 



(T) 



1 



4o3^2 /-oo 



72900 m2 

dAa 



dR 



{Aa{R)yMR) 



2 „ . / dAa 



fA{R) 



+ Aq- 



di? 



/5(i?)+Aa^/6(i?) + /7(i?) 



dR 

e-^^^'-^R\ 



(69) 



where the auxiliary functions fsiR) to fT{R) are defined by the following expressions: 



/3(^) 



480 104 fdV 



i?4 



176dyd2y 
~R~ciRchR2^ 

,dy dW 



44 



U) di^^' + ^^di^ 



96 d^V" 



/3 (70) 



+ 48 



di? d/?3 



^ dR^^ 
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3. Quantum-statistical theory of the Cotton-Mouton effect 



The birefringence of the refractive index can also be observed in the magnetic field B, and 
the relevant quantity describing this effect is the Cotton-Mouton constant Cm defined by 
Buckingham and Pople by the following expression [5]: 

^™=3(n2 + 2)2|^o~B^' ^^^^ 

In analogy to the Kerr constant K^, we can write the following virial expansion of the 
Cotton-Mouton constant: 



Cm = Ac MP + BcM (r)p2 + . . . , (76) 
where Aqm is proportional to the atomic electric-magnetic second hyperpolarizability r/g: 

AcM = ^^0, (77) 

and the second Cotton-Mouton virial coefficient is given by an expression analogical to Eq. 
(16): 



BcM = ^ 



d^B2{T-F\\-B) d^B2{T-F^-B) 



dB^dF^ dB^dF^ 



(78) 

B=F«=F^=0 



The two-body Hamiltonians describing the relative motion of two atoms in the mixed 
electric i^li/-*- and magnetic B fields are given by: 

^(F||,B) = m - kfzim' - ^«rz(R)i^if - ]4izzim'Fi (79) 



H{F^,B) = m - l^fAm' - lafA^)Fl - \v¥z,xxim'Fl (80) 

^(0) 

is given by Eq. (24). The space-fixed quantities ^, a, and rj are collision-induced mag- 
netizability, polar izability, and mixed electric-magnetic hyperpolarizability, respectively. 
They can be related to the body-fixed quantities by the expressions identical to Eqs. (27)- 
(28) and (32)-(33). Similarly as for the Kerr virial coefficient, it is useful to split the expres- 
sion for the Cotton-Mouton virial coefficient into contributions due to the electric-magnetic 
hyperpolarizability BcM,rj{T) and the anisotropies of the collision- induced polarizability 
and magnetizability BcM,/\£_Aa{T)'- 

Bcm{T) = BcMAi^a{T) + BcM,v{T). (81) 
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The quantum-statistical expression for BcM-q{T) reads: 

1,3 ^ / (_^\J+2I 



(82) 



Y,e-P^-^{i,n.j{R)HR)\^nj{R)) 
n 

+ / dEke-f'^^i^PE^jiRMRME^jiR)) 
■Jo 

where r]{R) is given by the following combination of the body-fixed Cartesian components: 



v{R) =— 

' 15 



'^VxXjXxiR) ^^xx,yy (R) - 2r],,^,,iR) 
+ l2r]„{R) - 27],,^^^{R) + 27],,^,,{R) 
The expression for BcM,Ai,AaiT) has the following form: 



(83) 
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+ {i;E,J'{R)\AC{RMnJ{R)){i^nJ{R)\^a{R)\^^;E,J'{R))) 



+ / dE. / dE,, 



El. — El. 



X {i;EM\^O^iR)\i^E,,J'iR)){i^E,,J'iR)\^CiRME.jiR)) 



(84) 



We end this section by saying that the semiclassical expansion for the Cotton-Mouton 
virial coefficient can be obtained in the very same way as described in sec. 2.4 for the Kerr 
virial coefficients, so we do not report explicit expressions here. 



4. Numerical results and discussion 



4.1. Computational details 



All numerical results reported in this section were obtained for the bosonic ^He isotope. 
The interatomic interaction potential was taken from Ref. [21], while the anisotropy of 
the collision- induced polarizability tensor from Ref. [14]. All remaining collision- induced 
properties were taken from the works of Rizzo and collaborators [22, 23]. It should be 
stressed here that strictly speaking the electric and mixed electric-magnetic properties 
should be taken at the frequency of the non-resonant fields and F-^. However, as shown 
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in Ref. [14] the frequency dependence of the polarizabihty is very weak in the frequency 
range used in the experiment, so the frequency dependence can safely be neglected. 

The Schrodinger equation for the relative motion was solved with the de Vogelaere 
method [32] which allows accurate calculations of the wave functions with an error quar- 
tic in the integration step, and computationally less demanding than other fourth-order 
algorithms, e.g. the Numerov method [33]. The matrix elements of the collision-induced 
properties with the radial wave functions were computed with the generalized Simpson 
method with the convergence criterion of a relative error of 10~^. The integration was 
done on the interval from 3 to 200 bohr. For numerical convenience, the integration over 
the energy Ej. was replaced by the integration over the wave vector k. Integration over k 
was done in the range 0.01 to 15 a.u. with a step of 0.002 a.u. We have checked that the 
contribution from the high k region, with k above 10 a.u., was very small. 

Numerical calculations of contributions to the Kerr virial and Cotton-Mouton virial 
coefficients related to the anisotropy of the collision-induced polarizabihty (and magneti- 
zability) are somewhat more complicated since they involve a double integration over k. 
The functions under the integral sign in Eqs. (56) and (84) have singularities at = Ek'- 
This singularity can be removed by using the following identity: 

lim = (3e-^'^', (85) 

but it may be the source of potential numerical inaccuracies. Therefore, it is advantageous 
to use the following integral representation of this singularity: 



Ek — Ek' 



e 



/ d(Te-'^(^^'-^'=). (86) 







The additional integration over a does not introduce any significant complications in our 
numerical procedure, as it does not require calculations of any additional matrix elements 
of the type (V'^fc j| Aa|'i/'_Bj^, j'). The latter calculations represent by far the most consuming 
step in our numerical procedure. The number of partial waves in the summations was such 
that the final result was converged within 1% at worst. 



4.2. Quantum- statistical results for the Kerr virial coefficient of the helium-^ gas 

We start the discussion with the analysis of the coUisional hyperpolarizability contribution 
to the second Kerr virial coefficient as a function of the temperature T. The results of the 
quantum statistical calculations of BK,-y{T) as function of the temperature are presented 
in Table 1 and illustrated in Fig. 2. Also presented in this Table is the classical term com- 
puted with Eq. (62), B^^^^(T), and the first and second quantum corrections, B^^^^{T) and 

B^^^{T), respectively. The quantum corrections have been computed from the expressions 
reported in sec. 2.3. An inspection of Table 1 shows that the quantum effects are small 
for temperatures larger than 100 K, and Bp^yiT) can be approximated by the classical 
expression with an error smaller than 5%. At lower temperatures the hyperpolarizability 
contribution to the Kerr virial coefficient of the '^He gas starts to deviate from the classical 
value. Still, for T > 50 K the quantum effects can efficiently be accounted for by the sum 
of the first and second quantum corrections. Indeed, for T = 50, 75, and 100 K the series 
Bf]^ = bP^^{T) + B^^^^^{T) + B^^^^^{T) reproduces the exact results with errors smaller than 
2%. One may note that at these temperatures the second quantum correction is small, and 
can be neglected for all practical purposes. In fact, the sum bP^{T) + bP^{T) slightly 
overestimates the exact result, while the full semiclassical term through the second order. 
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Table 1. Collisional hyperpolarizability contribution to the Kerr virial coefficient of helium-4 (in atomic units) as a function 
of the temperature T (in K). The consecutive columns report the temperature T, the classical result _B^^_^, the first quantum 
correction B^-*^, the second quantum correction scmiclassical result -B^^^, the Padc approximant [1/1], and the full 

quantum result BKn- 
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Figure 2. Second Kerr virial coefficient for the fielium-4 gas as a function of the temperature. Contribution i3x,7 
from the coUisional hyperpolarizability (in atomic units). 

7(^)1 slightly underestimates it. At temperatures below 50 K the semiclassical expan- 
sion in powers of starts to diverge. This divergence is clearly illustrated in Fig. 2, where 
the semiclassical result B^^^{T) and full quantum result BK,'y{T) are plotted as a function 
of the temperature. This behaviour of the power series in h? is not surprising, since the 
semiclassical expansion of the pressure virial and dielectric virial coefficients are known to 
diverge as well (see Refs. [9, 12]). 

Given the overall pattern of convergence of the semiclassical expansion, it is interesting 
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to find whether any rational approximations involving the low-order quantum corrections 
will reproduce the converged quantum result. It is well known [34-36] that divergent series 
can be effectively summed by means of Fade approximants. Since we know only three terms 
in the expansion of BK,y{T) as a power series in h'^, we could only use the simplest [1/1] 
approximant defined by the following expression: 



[1/1] 



?(0) 



>(2) 



B^'iT) - B\^>{T) 



(2) 



(87) 



The values of this approximant at various temperatures are reported in the sixth column 
of Table 1. Except for the lowest temperatures, the simple [1/1] Fade approximant works 
surprisingly well. For T = 15 and 20 K the sum of the classical term and first and second 
quantum corrections overestimates the exact result by 211% and 38%, respectively, while 
the [1/1] approximant reproduces the quantum results with errors of the order of 5%. This 
result is very gratifying since the calculation of the quantum corrections is much simpler 
than full quantum-statistical calculations. It is worth noting here that similar results were 
obtained for the second dielectric virial coefficient of the helium-4 gas [12]. 
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Figure 3. Second Kerr virial coefficient for the hclium-4 gas as a function of the temperature. Contribution Bk,a. 
from the polarizabihty anisotropy (in atomic units). 



We continue the discussion with the analysis of the contribution due to the anisotropy of 
the collision-induced polarizabihty tensor to the second Kerr virial coefficient as a function 
of the temperature T. The results of the quantum statistical calculations of B^/xaiT) 
as function of the temperature are presented in Table 2 and illustrated in Fig. 3. Also 
presented in this Table is the classical term B^^^^^{T), and the first and second quantum 

corrections, B^j^^^^{T) and -B^''^^(T), respectively. The quantum corrections have been 
computed from the expressions reported in sec. 2.3. An inspection of Table 2 shows that 
also in this case the quantum effects are small for temperatures larger than 100 K, and 
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Table 2. Collisional anisotropy contribution to the Kerr virial coefficient of helium-4 (in atomic units) as a function of the 
temperature T (in K). The consecutive columns report the temperature T, the classical result ^^^z^^, the first quantum 
correction B^j^^^^, the second quantum correction -B^^^^, scmiclassical result ^x^Aq' Fade approximant [1/1], the 
approximate result according to the prescription of Rizzo and collaborators [15] AcT ' ^^^'^ ^^'^^ quantum result B^./Sa ■ 
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BK,Aa{T) can safely be approximated by the classical expression with an error smaller 
than 2%. At lower temperatures the quantum result starts to deviate from the classical 
value. Still, for T > 50 K the quantum effects can very efficiently be accounted for by the 
sum of the first and second quantum corrections. Similarly as in the case of the collisional 
hyperpolarizability contribution, for T = 50, 75, and 100 K the series B^^^^ = B^^ ^^{T) + 

B^j^^^{T) + B^^^^[T) reproduces the exact results with errors smaller than 1%. Again, 
at these temperatures the second quantum correction is small, and can be neglected for 
all practical purposes. Around the temperature of 20 K the semiclassical expansion in 
powers of fi? starts to diverge. This divergence is clearly illustrated in Fig. 3, where the 
semiclassical result B^^^^{T) and full quantum result BK,Aa{T) are plotted as a function of 
the temperature. Similarly as in the case of i?K,7 the performace of the simplest [1/1] Fade 
approximant is very good. The divergent semiclassical series can effectively be summed up 
at temperatures as low as 15 K, and even at 10 K the error of the [1/1] approximant with 
respect to the exact quantum result Bx,AaiT) does not exceed 6%. Finally, we note that 
the approximate expression for BK^Aot{T) (denoted as i3^^°^(T)) advocated by Rizzo and 
collaborators [15] does not do a good job. Actually, at temperatures higher than 30 K the 
performance of this approximate quantum expression is worse than of the semiclassical 
expression. At lower temperatures, up to T = 10 K, the Fade approximant reproduces the 
full quantum result with a better accuracy than i?^^°^(T) 

Finally, in Table 3 and Fig. 4 we show the full Kerr virial coefficient, Bk{T) = 
BK,'yiT) + BK,Aa{T) as a function of the temperature, the classical term B^^^ , the semi- 
classical approximation and the approximate quantum result based on the expression 
of Rizzo and collaborators [15], B^^^°^. An inspection of Table 3 and Fig. 4 shows that 
the second Kerr virial coefficient is a smooth function of the temperature. It monotoni- 
cally decreases with the temperature, and around the room temperature it crosses zero 
and becomes negative. The overall performance of the classical approximation and of the 
semiclassical expansion are approximately the same as for the contributions Bx^-yiT) and 
BK,Aa{T). We note again that the semiclassical expansion gives more accurate results than 
the approximate expression of Ref. [15] for temperatures as low as 30 K. Only below 30 K 
^approx ]-,g(,Qj^gg closer to the exact quantum result Bk, but the error is large, 10%, 13%, 
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Tabic 3. Second Kerr virial coefficient (in atomic units) as a function of the temperature (in K) . Comparison of the quantum 
with the approximate results according to the prescription of Rizzo and collaborators [15] ^/^AcT' scmiclassical results. 
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and 17% for T = 20 K, 15 K, and 10 K. Thus, this expression does not seem to be useful, 
at least for the helium gas. 
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4.3. Comparison with the experimental data 



The Kerr constant Km as defined by Buckingham [2, 3], Eq. (3), is not measured experi- 
mentally. However, the result of the experimental measurements is: 



exp 



2 , n 
— lim — 
27 £^ 



n± 



AkP+ Bk{T) + Ak 



27ran 



.2 + 0(p3), 



where Ak is given by Eq. (19) and ao stands for the atomic polarizability. Eq. (88) shows 
that the coefficient multiplying p'^ is composed of two terms: the Kerr virial coefficient 
depending on the temperature T, and a T independent term. The latter can be evaluated 
by taking ao = 1.383181 a.u. and 70 = 43.1042 a.u. [37]. A simple arithmetics shows that 
even at high pressures of the order of 300 kPa the T independent term is very small, and 
the major part of the quadratic behavior must be related to Bk{T). 

However, despite a few attempts [6, 38, 39] the second Kerr virial coefficient for the 
helium gas could not be measured. In all experiments reported thus far, the dependence of 
the Kerr constant on the gas number density p was linear. Only the Authors of Ref. [39] 
tried to estimate the contribution from the pair interactions to the helium Kerr effect, 
however due to the limited pressure applied in the experiment the uncertainty of the fitted 
Bk at T = 294.8 K is huge, —366 it 499 a.u.^, making the result not very useful for 
comparisons between theory and experiment. All the measurements were carried out in 
the temperature range between 240 and 300 K, so in view of our results difficulties with 
observation of the two-body effects in the Kerr experiment are not surprising, since in this 
range of temperatures B^ is very small for helium and around 300 K it even crosses zero. 
A simple estimate based on our results demonstrates that for the gas number density of 
the order of 10^^ atoms per cm^ the contribution of BkP^ to K^y^p becomes important for 
temperatures below 100 K. Assuming the experimental precision of the order of 5%, actual 
observation would be possible at temperatures of the order of T = 10 K or below. 

In 2004 a measurement of the Kerr constant -ftTexp was reported for the superfluid helium 
in the temperature range 1.5-2.17 K [40]. The measured value of -RTexp = (1-43 it 0.06) x 
2Q-20 (|cni/V)^ can be compared with our calculations. The atomic contribution at the 
liquid helium density can be estimated to be A^p = 1.10 x lO"^*' (cm/V)^. The difference 
of 0.33x 10"^'^ (cm/V)^ is due to pair and nonadditive three-body (and higher) interactions. 
On the basis of our results in this temperature range the contribution from the second Kerr 
virial coefficient Bk{T) should be of the order of (5.18 — 3.88) x 10"^*^ (cm/V)^, which 
differs by an order of magnitude from the experimental result. However, such a difference 
between the liquid phase and gas phase results is not very surprising and was observed in 
many cases [41]. 

Some information on the Kerr virial coefficient for an atomic gas can be obtained from 
the analysis of the depolarized Raman spectrum. At high temperatures the lowest moment 
of this spectrum Mq defined as: 



Mn = - — 




(89) 



where /y is the intensity of the depolarized band as a function of the frequency u and Aq 
is the wavelength of the laser light, is proportional to the part of the Kerr virial coefficient 
related to the anisotropy of the interaction- induced polarizability tensor Bx,Aa{T) by the 



The applied conversion factor for Bj^ from atomic units to SI is 1 a.u.= 3.01154 10 C^m^J ^. 
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Figure 5. Total second Cotton-Mouton virial coefficient (in a.u.) for the hclium-4 gas as a function of the temperature 
(in K). 

expression [13]: 

Note that this relation is strictly valid only in the limit of high temperatures. At low 
temperatures it can only be related to the upper bound of Bruch et al. [13, 16]. However, 
at high temperatures we can compare the present theoretical value with the measured 

zeroth moment of Ref. [25]. The experimental value of Mq = (1.70±0.19) x 10"^ [25] is 
equivalent to Bx,Aa{T) = 86.9 it 9.7 a.u. which compares relatively well with the present 
theoretical result of 101.4 a.u. The measured value of Mq at low temperature of 99.6 K 
, Mq = 1.46 X 10~ A , can be translated to Bx^Aa = 220.7 a.u. computed from the 
expression of Bruch et al. [13, 16]. The computed value from this approximate expression 
is 265.7 a.u. It is noticeable that the computed values are systematically higher by ~ 20% 
from the results derived from the experiment. Given the fact the spectral moment Mq 
is obtained from the integration of the experimental intensity of the depolarized band, 
which is affected by some background intensity that has to be eliminated, such an 
agreement between theory and experiment should be considered as satisfactory. 



4.4. Cotton-Mouton effect in the helium-^ gas 

The second Cotton-Mouton virial coefficient as a function of the temperature T in the 
range from 4 K to 323 K is reported in Table 4 and graphically illustrated in Fig. 5. As 
inspection of Table 4 and Fig. 5 shows that the Cotton-Mouton virial coefficient is a smooth 
function of the temperature slowly increasing with T. At all temperatures considered in 
the present paper Bqm is negative. Also reported in Fig. 5 are the contributions to Bcm 
from BcM,AiAa{T) and BcM,rj{T)- Similarly as in the case of the Kerr virial coefficient 
at low temperatures the contribution from BcM,A^Aa{T) largely dominates. Only around 
200 K BcM,A^Aa{T) and BcM,n{T) become equal, and at room temperature BcM,-q{T) is 
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Tabic 4. Classical C^^^,^), semiclassical (B^j^j) and quantum ( Bcm) results for the second Cotton-Mouton virial coeffi- 
cient (in a.u.) for the helium-4 gas as functions of the temperature (in atomic units). 



T 


n(0) 


rjscl 


Bcm 


4 


-564.85 


-37925.98 


-43.22 


7 


-131.11 


-1388.97 


-32.03 


10 


-69.10 


-248.00 


-26.68 


15 


-40.12 


-56.02 


-21.97 


20 


-29.80 


-29.18 


-19.34 


30 


-21.53 


-18.33 


-16.47 


40 


-18.02 


-15.50 


-14.90 


50 


-16.07 


-14.16 


-13.92 


75 


-13.69 


-12.60 


-12.55 


100 


-12.60 


-11.88 


-11.87 


150 


-11.64 


-11.23 


-11.23 


200 


-11.27 


-10.99 


-10.99 


250 


-11.11 


-10.91 


-10.91 


300 


-11.07 


-10.91 


-10.91 


323 


-11.08 


-10.93 


-10.92 



larger, although not dominant. 

In Table 4 we also report the classical term and the sum of the semiclassical expansion 
through the second order. Similarly as in the case of the Kerr virial coefficient the classical 
expression works well up to the temperatures of wlOO K or higher. At lower temperatures 
quantum effects start to play the game, but still the semiclassical expansion effectively 
takes into account quantum effects at temperatures as low as 40 K. Indeed, at T = 40 K 
the error of the semiclassical result with respect to the full quantum result is only 4%. At 
still lower temperatures it diverges very fast. For instance, at T = 20 K the semiclassical 
result overestimates (in the absolute value) the quantum result by 51%, and at 15 K its 
absolute value is almost three times larger. 

There were a few measurements of the Cotton-Mouton effect in helium-4 gas [42-45] in 
the temperature range between 285 K and 300 K, i.e. close to the room temperature, and 
at relatively low gas number densities (the pressure of 1 atm, and gas number density of 
the order of 10~^*^ atoms per cm^). Measurements in Refs. [42, 43] were performed for a 
single gas pressure, assuming that a linear dependence of the Cotton-Mouton effect on the 
gas density is fulfilled. Also in these experiments the uncertainty of the measurements was 
relatively high, of the order of 20%, which precluded observation of any fine effects due to 
the interatomic interactions. The most accurate experimental data for the Cotton-Mouton 
effect in helium were reported in Ref. [44]. Although the precision of measurements was 
much higher than in the previous experiments it was only possible to observe the linear 
term in the virial expansion (76). On the basis of our results we predict a contribution 
of the Cotton-Mouton virial coefficient to the Cotton-Mouton constant Cm between 3.7 
ppm at the room temperature to 14.8 ppm at T = 4 K for the largest gas number density 
considered in Ref. [44]. Thus, it seems that similarly as in the case of the Kerr effect, 
experimental observation of the effect of interatomic interactions on the Cotton-Mouton 
constant for helium gas will be very challenging. 



5. Summary and conclusions 

The results reported in the present paper can be summarized as follows: 
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(1) Theory of the birefringence of the refractive index in atomic diamagnetic dilute gases 
in the presence of static electric (optical Kerr effect) and magnetic ( Cotton- Mouton 
effect) fields was formulated, and virial expansions of the Kerr and Cotton-Mouton 
constants as power series in the gas number density p were derived. It was shown that 
both virial coefficients can rigorously be related to the difference of the fourth deriva- 
tives of the thermodynamic (pressure) virial coefficient with respect to the strength of 
the non-resonant optical fields with parallel and perpendicular polarizations and with 
respect to the external static (electric or magnetic) field. Explicit quantum-statistical 
expressions for the second Kerr and Cotton-Mouton virial coefficients valid both in 
the low and high temperature regime in terms of the collision-induced electric and 
magnetic properties, and eigenvalues and eigenfunctions of the field free Hamilto- 
nian describing relative nuclear motion of two interacting atoms were derived. Our 
quantum-statistical expressions are significantly different from the approximate ex- 
pressions reported by Rizzo and collaborators [15]. 

(2) Semiclassical expansion of the second Kerr virial coefficient as a power series in fi? 
was derived, and explicit expressions for the first and second quantum corrections 
to the classical expression were reported. The consecutive terms in the semiclassi- 
cal expansion were expressed as one-dimensional integrals involving the Boltzmann 
factor exp(— /31/) and some functions depending on the potential, collision-induced 
properties, and their derivatives with respect to the interatomic distance R. 

(3) Both the second Kerr and Cotton-Mouton virial coefficients are smooth functions of 
the temperature. The Kerr virial coefficient monotonically decreases with the tem- 
perature and around the room temperature it crosses zero and becomes negative. The 
Cotton-Mouton virial coefficient is also a monotonic function of the temperature. In 
the range of temperatures considered in the present paper it is always negative and 
slowly increases with T. 

(4) Semiclassical expansions through the second order are shown to diverge at low tem- 
peratures. However, for helium gas they account for all major quantum effects up to 
the temperatures of ~50 K. At T = 50 K the errors of the semiclassical expansion of 
the Kerr and Cotton-Mouton virial coefficients with respect to the quantum results 
are 0.7% and 1.7%, respectively. The simplest Fade approximant greatly improves 
the convergence, and effectively sums up the divergent series at temperatures as low 
as 20 K. In the temperature range where the semiclassical expansion is valid, the 
semiclassical and full quantum results agree better, than the approximate results 
according to Ref. [15] and present quantum results. 

(5) Despite many efforts, neither the Kerr nor Cotton-Mouton virial coefficients could 
be measured for the helium gas thus far. Based on the present numerical results, 
estimates of the temperature range for which the effects could be observed were 
reported. They mostly concern very low temperatures that are not easily accessible 
to the gas phase experiments. It seems that experimental data for other light atomic 
gas like neon are necessary to judge the importance of the quantum effects on the 
optical birefringence at very low temperatures. 
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